DECOMPOSITION OF WAVELET REPRESENTATIONS AND MARTIN 

BOUNDARIES 
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Abstract. We study a decomposition problem for a class of unitary representations associated 
with wavelet analysis, wavelet representations, but our framework is wider and has applications to 
multi-scale expansions arising in dynamical systems theory for non-invertible endomorphisms. 

Our main results offer a direct integral decomposition for the general wavelet representation, 
and we solve a question posed by Judith Packer. This entails a direct integral decomposition of the 
general wavelet representation. We further give a detailed analysis of the measures contributing to 
the decomposition into irreducible representations. We prove results for associated Martin bound- 
aries, relevant for the understanding of wavelet filters and induced random walks, as well as classes 
of harmonic functions. 

Our setting entails representations built from certain finite-to-one endomorphisms r in compact 
metric spaces X, and we study their dilations to automorphisms in induced solenoids. Our wavelet 
representations are covariant systems formed from the dilated automorphisms. They depend on 
assigned measures /i on X. It is known that when the data {X, r, jj.) are given the associated wavelet 
representation is typically reducible. By introducing wavelet filters associated to {X, r) we build 
random walks in X, path-space measures, harmonic functions, and an associated Martin boundary. 

We construct measures on the solenoid {Xoo,roo), built from {X,r). We show that Voo induces 
unitary operators U on Hilbert space H and representations tt of the algebra L°°{X) such that 
the pair (?7, roo), together with the corresponding representation n forms a crossed- product in the 
sense of C*-algebras. We note that the traditional wavelet representations fall within this wider 
framework of ("H, U, tt) covariant crossed products. 
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1. Introduction 

We study a decomposition problem for a class of unitary representations associated with wavelet 
analysis, even though our framework is wider and has applications outside multi-scale wavelet 
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expansions, see details below. One powerful tool in the construction of families of multi-scale 
wavelets (see e.g., |Dau921 [Mal98l [Jor06] ) is an introduction of a finite system of filters. Here we 
understand the notion of "filter" in the sense of signal processing. In this context, each filter will 
be a function of a complex number z (a frequency variable), and for many purposes it is enough 
to consider only a phase of 2, so we may restrict attention to the case when points z are in the 
1-torus. Each function mi{z), i = 0,1, . . . , N — 1 typically supports a frequency band. With these 
conventions, in wavelet considerations, the function mo represents a low-pass filter, i.e., passing 
low frequency signals. When a suitable Fourier expansion is introduced for the filter functions we 
arrive at the masking coefficients that determine some particular wavelet. This framework includes 
both traditional wavelet systems in the Hilbert space L^(M'^) in some number of dimensions d, 
as well as orthonormal wavelet bases on fractals, as studied by two of the present authors, see 

[Djn7llDjn6liBKinj . 

Analysis of filters. Continuing with filter functions tuq on the circle group T, we consider, for 
every fixed z in T, then the absolute squared mo with some normalization W{z) = \mQ{z)\'^/N . 
With this we then get a family of probability distributions on the set of N solutions in T to the 
equation = z (see p.ip - p.5p below). There will be a solution w in each of the A'^ frequency- 
bands. In the special case when N = 2, the function mo, or the system of functions mo and 
mi, are called a quadrature mirror filter (QMF). The reason for this is that two other operations, 
down-sampling, and up-sampling, allow one to build a discrete wavelet algorithm with dual filters, 
the dual one is the "mirror". 

Here we have adopted a more general framework: Instead of T we will consider a compact metric 
space X, and our filter functions will be functions from X into the complex plane. We further 
generalize the choice of endomorphism r{z) = z^ , considering here instead an endomorphism r in 
X which is onto, and for which each pre-image of points in X is a finite subset; so finite-to-one 
endomorphisms. In both the traditional wavelet case, and in the more general framework, we end 
up with dynamical systems in a solenoid. 

When the pair {X, r) is given as specified, there is a standard way of building a solenoid Xqo = 
-^00 ('^) over X (see ()2.7p and (j2.8p below). There are several advantages working with the solenoid: 

(i) The endomorphism r in X induces an automorphism r^o in Xoo(see (12. Sp ). 

(ii) The transition probabilities W on X induce a random walk on the solenoid which encodes 
properties of the representations induced by the prescribed wavelet filters [DJ07j : these 
representations are known as wavelet representations in the literature. 

(iii) With this random walk we are able to compute transition probabilities, harmonic functions, 
and associated Martin boundaries. 

(iv) Points in the solenoid may be thought of as random-walk paths; for each point x in X, we 
will have an infinite random-walk path, represented as a subset of the solenoid. 

(v) Fixing the function W = \mQ{-)\'^ /^r^^{r[-)) on X we get a path space measure P. We 
will be interested in the family of measures {Px), x m. X, with P^ conditioned on the set 
of paths starting at x (see (j2.12p ). 

To specify wavelet representations we must then also have a prescribed measure // on X (see 
Definition 12. ip . For the theory of wavelet representations, see for example [HL081 ILar07bl ILarOTaj 
IPac041 IPacOSbl IPacOSa] . For an early treatment of transfer operators Rw in multi-scale wave- 
lets, see |Jor01| . The relevance of the Rw harmonic functions (i.e., solutions h to Rwh = h) in 
the decomposition theory for the corresponding wavelet representation p\Y was pointed out there. 
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Specifically, [JorOlj has the idea of computing operators in the commutant of pw from i?pi/-harmonic 
functions. 

For computations of i?vi/-harmonic functions of concrete wavelets, see also |BJ02j . The theory 
of Rw harmonic functions was developed recently by a number of other authors; in |Dut06] is 
introduced the study of intertwining operators for pairs of wavelet representations. 

In our present wider context, the i?iy-harmonic functions enter in equation (j2.5p below, and they 
underlie our considerations in Section 4 regarding Martin boundary. 

In our present setting, in an earlier paper [DS] a necessary condition was given on the data 
{X, r, n) for when the associated wavelet representation is reducible. For the reader's convenience, 
we have stated it as Theorem 12.41 below. 

Summary of results. Our main result here (Theorem 13. 3p offers a direct integral decomposition 
for the general wavelet representation. This completely solves a question posed by Professor Judith 
Packer, see e.g., [PacOSal iBFMPlnl iBLP+lol iBLM+OOl IBFMPOQj . In our Theorem 13.31 we offer 
a direct integral decomposition of the general wavelet representation, and our Theorem 14.51 deals 
with a derivation of the measures contributing to the decomposition. Our results yield as clean a 
decomposition for the general wavelet representation into irreducibles as is realistically feasible. 

In Section 4, we have included a result on an associated Martin boundary. Even though it is not 
used directly, it is certainly relevant for the understanding of our random-walk harmonic functions, 
and the question about the Martin boundary naturally presents itself. In fact, we obtained our 
result in Section 4 in response to a question asked of us by Erin Pearse. 

Indeed, there are some intriguing connections to random- walk models studied recently in papers 
by one of the present authors and E. Pearse (see e.g., [JPlOj and [DJlOj . These are computations 
for infinite weighted graphs G. In both instances we get transition probabilities and associated 
random walks. 

If G is a graph as in [jPToj, the condit ion for the context of these studies is that, for every vertex 
in G, there are only a finite number of transitions possible to neighboring vertices. But the Markov 
processes in |JP10] are reversible, and therefore the associated boundaries are more amenable. By 
contrast, our transition processes are non-reversible, except for some trivial special cases. More 
specifically, the reason our random walks are not reversible, is that transition happens from one 
point X to one of the distinct subnodes y (neighbors) were y will be one of the finite number of 
solutions to the equation r(y) = x. Hence these transitions never return, unless we consider cases 
when X might be a periodic point, in which case they might return. 

2. Measures on the solenoid 

Before turning to our direct integral result, we begin with some preliminaries regarding measures 
on solenoids. Since our starting point is a given finite-to-one endomorphism r in a compact metric 
space X, it is then natural to look for a way of corresponding to this a unitary operator [/ in a 
Hilbert space H, such that U together with {X, r) satisfy a covariance relation; see (i) in Theorem 
12.21 below. The introduction of suitable measures on the associated solenoid (Xoc^oo), built from 
(X, r), then gets us a representation vr of the algebra L°°(X) such that [/, together with Too, form 
a crossed-product in the sense of C*-algebras. This is possible since is an automorphism. We 
will refer to a crossed-product system {T-L, UjTt) as a wavelet representation. 

Indeed, in |DJ07j . we proved that the traditional wavelet representations fall within this wider 
framework of {T-L,U,7r) covariant crossed products. Specifically, in the special case when X = T, 
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and the endomorphism r is just the power mapping r{z) = (for a fixed integer N > 1), then 
it can be seen that a covariant crossed products indeed speciahzes to a unitary representation of 
a corresponding A^-Baumslag-Sohtar group; see e.g., |DJ081 IDut06j . Even in the case of these 
classical Baumslag-Solitar groups, our understanding of the unitary representations and their de- 
compositions is so far only partial. 

Definition 2.1. Let X be a compact metric space and r : X ^ X he a finite-to-one, onto, Borel 
measurable map. Let /i be a strongly invariant Borel probability measure on X, i.e. 

(2-1) 1^"^^ = ! #^ ^ fiy)M^), 

r{y)=x 

for any bounded Borel function on X. 

A function tuq on X is called a quadrature mirror filter ( QMF) if 

(2-2) E K(y)P = l, (x€X) 

r{y)=x 

In what follows we will assume that: 

(2.3) the set of zeroes for tuq has /^-measure zero. 
Given a QMF vtiq we define 

(2.4) W{x) = J^^^'^f (xeX). 

ifr ^[r[x)) 

Then the function W satisfies the following equation: 

(2.5) ^(2/) = !' {xeX). 

r{y)=x 



Equation (|2.5p can be interpreted as an assignment of transition probabilities: the probability of 
transition from x to y G r~^(x) is equal to W{y). 
A function /i on X is called Ryy -harmonic if 

(2.6) Yl W{y)h{y) = h{x), {x e X). 

r{y)=x 

Theorem 2.2. [ DJ07] There exists a Hilhert space T-L, a unitary operator U onH, a representation 
vr of L°°{X) on H and an element cpofTi such that 

(i) (Covariance) UTi{f)U-^ = 7r(/ o r) for all f G L°°(X). 

(ii) (Scaling equation) Uip = 7r{mQ)Lp 

(in) (Orthogonality) (7r(/)(/?, ip) = J f for all f G L°°(X). 

(iv) (Density) {[/-"7r(/)99 [ n G N, / G L°°{X)} is dense in Ti. 

Moreover they are unique up to isomorphism. 

Definition 2.3. We call the system {7i,U,Tr,varphi) in Theorem 12.51 the wavelet representation 
associated to the function rriQ. 



WAVELET REPRESENTATIONS AND MARTIN BONDARIES 5 

We recall some facts from [DJ07| . The wavelet representation can be realized on a solenoid as 
follows. Let 

(2.7) Xoo := {(xo,xi,...) e X^|r(x 

n+i) — Xn for all n > oj- . 

We call the solenoid associated to the map r. 

On Xoo consider the cr-algebra generated by cylinder sets. Let r^o '■ X^o 

(2.8) roo(xo,2;i, . . . ) = (r(2;o), xq, xi, . . . ) for ah (xo,xi,...) G X^. 

Then is a measurable automorphism on X^o- 
Define 6*0 : X^o X, 

(2.9) 6'o(xo,xi, . . . ) = xq. 

The following commutative diagram summarizes the relation between the maps r,roo,OQ: 

X(yQ — ^ Xqq 

X ^ X 

Define for m > the projection 6m ■ X^o X, 

Gmip'^Oi -^1 ; • • • ) — •^m • 

The measure fioo on X^o will be defined by constructing some path measures Px on the fibers 

■= {ixo,Xi, . . . ) G Xoo I Xq = x}. 

On i^xo we will consider the infinite product topology which is defined by the basis of open sets: 
for n > 0, xi, . . . ,Xn & X with r(xj+i) = Xj, j G {0, . . . , n — 1}, 

(2.10) ^a;o,...,x„ := {(zo^Zi,...) G ^Ixo ■ Zq = Xq, . . . , Zn = Xn}- 

With this topology flxo is a compact Hausdorff space. 
Let 

c(x) := #r~^(r(x)), W{x) = \mo{x)\'^ /c{x), (x G X). 

Then 

(2.11) J2 W{y) = l, (xGX) 

r{y)=x 

W{y) can be thought of as the trasition probability from x = r{y) to one of its roots y. 
For X G X, the path measure Px on is defined on cylinder sets by 

(2.12) Px{{{Xn)n>0 enx\xi= Zi,...,Xn = Zn}) = W{zi) . . . W{Zn) 

for any zi, . . . ,Zn G X. 

This value can be interpreted as the probability of the random walk to go from x to Zn through 
the points xi, . . . , x^. 

Next, define the measure //qo on X^o by 

(2.13) fdfloo= / /(x,xi, . . . ) dPa;(x,xi, . . . ) d^(x) 

J Jx Jn^ 

for bounded measurable functions on X, 



oo ■ 
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Consider now the Hilbert space H. = Lp'{ii^). Define the operator 

(2.14) Ui = {mooeQ)ior^, e L\X^, fi^)). 
Define the representation of L°°(X) on T-L 

(2.15) 7rim = {fo0o)t if GL'^{X),^eL\X^,f,^)). 
Let (/9 = 1 be the constant function 1 on X^o. 

Theorem 2.4. [DJ07] Suppose tuq is non-singular, i.e., fj,{{x € X \mQ{x) = 0}) = 0. Then the 
data (1-1,17, Tr,ip) forms the wavelet representation associated to niQ. 

Theorem 2.5. jPS] Suppose r : {X,fj,) — > {X,n) is ergodic. Assume \mo\ is not constant 1 ^-a.e., 
non-singular, i.e., /i(mo(a;) = 0) = 0, and log |moP is in L^{X). Then the wavelet representation 
{T-L, U, vr, if) is reducible. 

We will be interested in the decomposition of the wavelet representation into irreducibles. We 
need a few more notations and lemmas. 

Definition 2.6. Define 

mo = 1, mn = (mo o 6q) ■ {mo o Oq o r^o) ■ ■ ■ {mo o 0o o r^^), for n > 1, 

rhn = rr 1 for n < 0. 

(mo o 00 o Too ) . . . (mo o 6*0 o r^,) 

The function m : X^q x Z —?■ C* defined by m{x,n) = mn{x) gives a one-cocycle for the action 
of Z on Xoo determined by Too- 

The fact that U is an isometry implies the following lemma. 
Lemma 2.7. For ^ E L'^{Xoo, fJ-oo): 

j (.dfioo = j Irhnl'^^, or^dfioo, (n G Z). 

3. The decomposition of the wavelet representation 

Definition 3.1. We say that a subset J- of X^o is a fundamental domain if, up to /ioo-measure 
zero: 

y r^(J-) = Xoo and r^{J^) D rZ{J') = forn / m. 

Definition 3.2. For z = {zo, zi, . . . ) in X^o define the following representation: consider the Hilbert 
space 



I neZ 



|?7T-„(z)P < OO 



with inner product 

(C, V)h. ■=^UVn\^n{z)\'^. 

Note that we avoid here the points z G X^a such that one of the functions mn{z) = 0. Since mo 
is non-singular, such points form a set of ^oo-measure zero. 
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Define the unitary operator 



Uz{Cn)nez = {mo o 6*0 o rj^(z)^„+i)„gz- 



Define the representation of n of L°°{X): 



'^z{f){in)n& = {fo6oo (z)e„)„6Z, (/ £ 



The representation vr^ is defined for bounded functions on X, not just essentially bounded. The 
//-measure zero sets will affect the individual representations vr^ but not their direct integral (see 
below) . 

Theorem 3.3. In the hypotheses of Theorem \2.5l there exist a fundamental domain T . The wavelet 
representation associated to mo has the following direct integral decomposition: 



where the component representations {Hz^Uz^t^z] in the decomposition are irreducible for a.e., z in 
T , relative to fioo- 

Proof. We state the irreducibility of the component representations in a lemma: 

Lemma 3.4. For almost every z € X^, the objects [Hz^Uz^'Kz] form an irreducible represen- 
tation. 

Proof. One has to check that Uz is unitary, vr^ is a representation and UzTTz{f)U~^ = "^zif ° ^) for 
all / € L°°{X). All these follow from simple computations. 

To see that the representation is irreducible for //oo-a.e. z, take z to be non-periodic, i.e., 
r^{z) ^ z for all n 7^ 0. Then {7r^(/) : / G L^{X)} forms a maximal abelian subalgebra with 
cyclic vector 5o (see |Tak021 Corollary III. 1.3]), where 5o{n) = 1 for n = 0, and 5o{n) = otherwise. 
Then, an operator A that commutes with Uz and tTz has to be of the form TTzig) for some g € L°°{X). 
Since A commutes with Uz we have iTzigor) = Uz'^z{9)U^^ = TTz{g)- This implies that g is constant 
on {r'^{z) : n G Z}, so A is a multiple of the identity. □ 

We begin the proof as in the proof of the main result in |DS| . 
From the QMF relation and the strong invariance of /i we have 



Since log is strictly concave, and |?tt,oP is not constant /u-a.e., it follows that the inequality is strict, 
and a < 0. 

Since r is ergodic, applying Birkoff 's ergodic theorem, we obtain that 





By Jensen's inequality we have 





n— 1 



k=0 



8 DORIN ERVIN DUTKAY, PALLE E.T. JORGENSEN, AND SERGEI SILVESTROV 

This implies that 

hm (jmo(x)mo(r(j;)) . . . mo(r"'~^(x))P)^''"' = e'^ < 1, /i — a.e. 

Take b with e" < 6 < 1. 

By Egorov's theorem, there exists a measurable set Aq, with )u(Ao) > 0, such that 

(Imo(x)mo(r(x)) . . . mo(r"~H^))P)^^" 

converges uniformly to on Aq. (Taking smaller if needed we can assume /^(^o) < 1-) This 
implies that there exists an uq such that for all m > uq: 

(|mo(x)r?T,o(r(x)) . . . mo(r"'"~^(x))p)^^™' < b for x € Aq 

and so 

(3.1) |mo(x)mo(r(2;)) . . . mo(r'"~^(x))p < 6™, for m > uq and all x G Aq. 

Next, given m G N, we compute the probability of a sequence {zn)nef>s G -'^oo to have Zm € ^o- 
We have, using the strong invariance of fi: 



P{Zm G ^o) = ^oo {{{Zn)n \ Zm G ^o}) = / 



' 41 -m( \ XI \mQ{zi)\^ ...\mQ{Zm)?XAo{Zm)d^l{zQ) 

r{zi)=zo,...,r{z,n}=z^^i 

\mQ{zm)mQ{r{zrn)) ■ ■ ■ mQ{r"'~'^{zm))fxAo{zm) dn{zm) 

X 

= / \mQ{x)mQ{r{x)) . . .mo{r"'~'^{x))\'^XAo{x) dfj.{x). 
Jx 

oo „ 

^ P(zm € vlo) = ^ / jmo(x)mo(r(x)) . . .mo(r"~^(x))pXAo < OO 

m=l m>l '''^ 



Then 



and we used p.ip in the last inequality. 

Now we can use Borel-Cantelli's lemma, to conclude that the probability that Zm G Aq infinitely 
often is zero. Thus, for /ioo-a-e. z := {zn)n, there exists kz (depending on the point) such that 
Zn ^0 for n > kz- In other words, if i?o = ^ \ ^0 then for /^oo-a.e. {zn)n in there exists k^ 
such that Zn € Bq for all n > kz- 

Define now the set 

^oo := {{zq, zi,...) e Xoo : ZQ,zi,--- E Bq} . 
It is clear that if {zq, zi, . . .) G Aoo then r^{zQ, zi, . . .) = (zi, Z2, . . .) is in Aqo too. Therefore 
r^{Aoo) ^ ^oo- This means also that A^o C roo(^oo)- 

From the statements above we see that for /Uco-a.e. {zq,zi, . . . ) in X^o there exists n such that 
Zn, Zn+i, ■ ■ ■ are in Bq which means that r^"'(zo, zi, . . .) is in A^o and so {zq, zi, . . . ) € r'^{Aoo)- 
Thus, up to measure zero: 

U r^(^oo) = Xoo. 
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We claim that also, up to measure zero, one has 

n c(^oo) = 0. 

Suppose {zq, zi, . . .) is in all r^{Aca) for n > 0. Then {r"'{zQ),r^~^ (zq), . . .) € Aqo so r^{zo) G 
Bq for all n > 0. Since < /i(i?o) < 1 this contradicts the fact that r is ergodic on X. 

Now take T := roo(^oo) \ ^oo- The properties of Aqo easily imply that is a fundamental 
domain. 

Next we check the direct integral decomposition. 
Define ^ : L'^iX^o, ^^oc) ^ f^n^ dfioo{z), 

We check that ^ is an isometry. We use Lemma 12.71 

= ljm){z)fnjf^oo{z) = \\nf- 

To check that ^ is onto, we can compute the inverse {^~^{S,{-)n)ne:z){z) = in{f^z) if z G rJ^(J'). 
Some direct computations show that ^ intertwines the [/-operators and the representations vr. 

□ 

4. Martin boundary 

The idea of associating to wavelet constructions a transfer operator Rw and associated harmonic 
functions was pioneered by W. Lawton in the two papers |Law90t ILaw91j . 

The idea is that wavelets are determined by a system of numbers, often called masking coeffi- 
cients. It is possible to turn these into coefficients in filter functions mj, and by selecting i = (see 
eq (12. 5p ) we get transition probabilities and a transfer operator Rw-, 

Rwf{x)= Yl i^^x). 

r{y)=x 

Hence Rw is determined by the prescribed masking coefficients, and the question is how prop- 
erties of the masking coefficients (and therefore of Rw) determine the wavelets. It turns out that 
this is decided by the spectrum of Rw, including the eigenspace for eigenvalue 1, which produces 
the harmonic functions. 

As shown in |DLS| . operators in the commutant of the wavelet representation correspond to 
bounded i?vi/-harmonic functions. If we restrict such harmonic functions to inverse orbits of points 
we get harmonic functions for the random walk, or what we call below p-harmonic functions. The 
Martin boundary theory offers a way to construct such harmonic functions by means of integrals 
on a certain boundary. We perform these computations here to see what the p-harmonic functions 
are in this case. 

Definition 4.1. A point xq € X is called regular if the following two conditions are satisfied: 

(i) The sets r~'^(xo), n G N are mutually disjoint. 

(ii) None of the sets r~"'(xo), n > intersect the set of zeroes of W. 
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Note that condition (i) means that xq is not periodic for the map r, i.e., r"'(xo) ^ xq for any n > 1. 

For a point xq E X, define the set T{xq) := U„>or~"'(xo). We cah this the tree with root at 
Xq. If Xq is regular and x G T(xo), define n{xQ) to be the unique non-negative integer such that 
r"(^o)(x) = Xq- 

Let Xq € X be regular. We define now a random walk on the set T{xq) and we construct its 
Martin boundary by following |Saw97j . 

For x,y T{xq) define the transition probabilities p{x,y) as follows: 

A function u on T{xq) is called p-harmonic if 

(4.2) u{x) = ^ p{x,y)u{y), (x G T(xo)). 

yeT(a;o) 

The function Pn{x,y) is the n-th matrix power oi p{x,y) and represents the probability of tran- 
sition from X to y in n steps. It can be easily seen that po{x, y) = 6xy and 

W{y)W{r{y)) . . . VF(r"-i(y)), if r"(y) = x, 



(4.3) Pnix,y) 



0, otherwise. 



The Green function or potential function is defined by 

oo 

(4.4) g{x,y) ■.= ^pn{x,y), (x,yGT(xo)). 

n=0 

Note that, for our random walk, only one term in the sum in (14. 4p is non-zero. 
The Martin kernel is defined by 

(4.5) Kix,y) := ^py^ (x,yGr(xo)). 

The denominator in ()4.5p is non-zero because each vertex y can be reached from xq eventually. 
The function K{x, •) is bounded by some constant Cx (which we will describe below). Set 

(4.6) M..y)= Y. (^^.y^TM). 

where D{q) > for all q G T(xo) and J2qeT{xo) ^ Here Sxy = 1 if x = y and = 
otherwise. 

The Martin compactification [7~(xo), p] is the completion of ^(xo) with the metric p. The Martin 
boundary is defined as dT{xo) := T(xo) \ T(xo). 

As shown in |Saw97j a sequence {yn} in T(xo) is Cauchy with respect to the metric p if and only if 
either (i) yn = U for all n > no for some y G T{xq) and some no G N, or else (ii) lim„_j.oo Vn = oo and 
linin^Qo K{x,yn) exists for all x G T(xo). (Here lim„_^oo = oo means that y„ leaves eventually 
any finite set and never returns.) 

Thus the Martin boundary 5T(xo) is the set of equivalence classes of Cauchy sequences that 
satisfy the condition (ii) above. 

The maps K{x,-), x G T(xo) extend uniquely to continuous maps on T(xo) and we use the same 
notation K(x, •) for their extensions. 
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Theorem 4.2. [Martin representation theorem] For any p-harmonic function u{x) > there 
exists a measure v on dTix^) such that 



(4.7) u{x) = / K{x,a)di'{a), (x G T(xo)). 

Jdnxo) 

Proposition 4.3. With the definitions above we have: 
(i) The Green function satisfies the equation 

W{y)W{r{y)) . . . W{r"-~^{y)), if r^{y) = x for some n > 0, 
0, otherwise. 



(4.8) g{x,y) 



(If n = the product is defined to be 1.) 

(ii) The Martin kernel is 

I XT.^/ / ^^ ^Tjr, r,t^\ 1 r ^^ =: C'z : if r^ (v) = X for somc n > 0, 

(4.9) K(x v) = i W{x)W{r{x))...W{r'^'-'^>-^{x)) x, j \yj j _ , 

' 1 0, otherwise. 

Thus K{x, •) is constant Cx on the subtree T{x) with root at x and everywhere else. Using 
the notation of Definition \2.6\ we have that K{x,y) = -^^^j^, if f^iv) = x; if #r~^{x) is 

constant, then the two functions Wn and rhn differ by a multiplicative constant (see (12. 4p ). 

(iii) A function u on T{xo) is p-harmonic if and only if 

(4.10) n(x)= ^ WivMy), (x G r(xo)). 

r{y)=x 

Proof, (i) follows directly from (14. 3p . Note that, because xq is regular, the number n such that 
^"(y) = a; is unique. For (ii), if r"'(y) ^ x for all n then g{x,y) = so K{x,y) = 0. If r"(y) = x, 
then we have 

g{x,y) W{y)...Wir'^-\y)) 1 



K{x,y) 



g{xo,y) iy(y)...VF(r"(y)-i(y)) W{r^iy)) . . .W{r^(y)-^{y)) 

1 



M^(x)... VF(r"W-i(x))' 
(iii) is obtained from the following computation: 

u{x) = ^p{x,y)u{y) = ^ p{x,y)u{y) = ^ W{y)u{y). 

y r{y)=x r{y)=x 

Theorem 4.4. Let xq be a regular point in X. Define the map $ : ~^ dT{xo) by 

(4.11) $(xo,xi,...) := {x„}, 

i.e., to each sequence in Q^o '^^ associate the equivalence class of this sequence in (?T(xo). 
Then $ is a bijective homeomorphism from ^Ix^ onto 9T(xo). 
For X G ^(xo) and (xo,xi, . . . ) G Qxo 

f w7 ^ \ — w7 7: if X = Xn for SOmC U > 0, 

(4.12) K(x,^(xo,xi,...)) = < W{xi)w{x2)...wix,,)^ J "-^ , 

[ 0, otherwise. 



□ 
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Proof. First, we show that $ is well defined, so is a sequence with the property that lim x„ = oo 
and lim K{x,Xn) exists for all x G T{xo). 

Since xq is regular, the sets r"^{xo) are disjoint, therefore any finite subset of T(xo) lies in a finite 
union Uj<jr~-'(xo), and since x„ E r~"(xo) for all n, it follows that x„ eventually leaves this finite 
set and never returns. 

For the second condition, take x € T{xq). Recall that n{x) is the unique number such that 
X G r~"(^^(a;o). We have two possibilities: x^^x) = x ov not. In the first case we have that x„ is 
in the subtree T{x) for all n > n{x) so K{x,Xn) is constant C^, by Proposition 14.31 In the second 
case, we have that Xn is not in the subtree T(x) for all n > n{x), so K{x,Xn) is constant 0. In 
both cases limi^(x, x^) exists, so $ is well defined. 

Next, we check that $ is onto. Take a sequence {yn} in Tixo) with limy„ = oo and such that 
limK{x, Un) exists for all x G T{xo). Since K{x, y) is either > or depending on whether y is 
in the subtree T{x) or not, it follows that for all x G T(xo), either the sequence is y„ is eventually 
contained in the subtree T{x) or it is eventually contained in the complement of T{x); it cannot 
jump back and forth between T{x) and its complement. 

We will construct by induction the sequence (a;„)„ in ilajg with ^{xq,xi, ■ ■ ■) = {yn}- The first 
element xq is given. Next consider the points zi,. . . ,zj in r~^{xo). Since limy„ = oo, eventually 
the sequence will be in the union of the subtrees L)j^^T{zj). With the previous remark, one of the 
sets T{zj) will contain the entire sequence yn eventually. We define xi to be the point Zj with this 
property. 

Inductively, if Xm has been defined such that the entire sequence y„ is in the subtree T{xm) 
eventually, we take the points in r~^{xni); since limy„ = oo, the entire sequence will lie in 
^zer-^{xm)'T{z) eventually. Since yn cannot jump back and forth between a subtree and its com- 
plement, there is one of the elements z G r~^{xni) such that the sequence y„ lies in the subtree 
T{z) eventually. We call this point x^+i- 

To prove that $(xo,xi . . . ) = {yn} we just have to show that the sequences {xn} and {yn} are 
equivalent, i.e., lim p{xn,yn) = 0. Take e > 0. There exists a finite subset F of T{xq) such that 
Ylq^F < e. For each q ^ F, either y„ is in the subtree T{yn) eventually, or it is in the 

complement of T{yn) eventually. From the definition of {xn} we see that x„ will have exactly the 
same property. Thus K{q, Xn) = K{q, yn) for n large enough. Also, since x„ and yn go to infinity, 
it follows that 6gx„ = Sgy^ for n large enough. Therefore, for n large, the terms in the sum in (14. 6p 
for p{xn,yn) that correspond to g G F are all zero, the rest are bounded by 2^g^^L'(g) < e. So 
pixn, yn) < e for n large, and therefore ^{xq, xi, . . .) = {yn} and <I> is onto. 

To see that $ is one-to-one, take / (x^) in ^Ixq- Let uq > 1 such that x„o ^ x'^^. Then for 
n > no, K{xno,Xn) = Cx and K{xno,x'n) = 0. Therefore 



p{xn,x'J > 



\-^i.XnQ ) Xf^ 



K{XnQ, X^)\ 



Cx„.. + 1 



a„„ + 1' 



This implies that ^{xn) ^ ^{x'n) so ^> is one-to-one. 

To prove that $ is continuous, take (x„) G rj^o and e > 0. Take a finite subset F of T(xo) such 
that '^'^q^p D{q) < e. Take no such that F is contained in Un<no'"~"(2;o)- Take (x^) G Kco,...,!:™^ 
so Xg = xo, • • . ,x^jj = x„Q. Then the sequences Xn and x^ are in the subtree Ti^Xn^) for n > no. 
This implies that for q £ F and n > no, we have that either both x„ and x^ lies in the subtree 
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T{q) or they both he outside T{q); also 6ga:„ = 6g^'^ = 0. Then < Y^g^p'^^il) < e so 

p{^{xn),^{x'n)) < e. This shows that <I> is continuous. 

Since both spaces Qxq and 5T(xo) are compact Hausdorff, it follows that <I> is a homeomorphism. 

Next, we check (j4.12p . Take x G T(xo) and (xq, xi, . . .) G Qxo- We have K{x, $(xo, xi, . . . )) = 
limi^(x, x„). If X ^ x„ for all n > 0, then Xn is not in the subtree T(x) for all n, so K{x,Xn) = 
for all n and therefore K{x, <J?(xo, xi, . . . )) =0. 

If X = x„ for some n, then for m > n, Xm is in the subtree T(x) so 

K{X,X„,) = K{Xn,Xm) = ^(^^)^(^(^^)) . . . ^(^n(.„)-l (^^)) = W {Xn)W {Xn-l) . . .W {x^)' 

This proves ()i32D . □ 
Theorem 4.5. For any p-harmonic function u > 0, there exists a measure v on fi^o such that 

(4.13) U{X) = ^(^)^(^(^)) ^ Vt^(^nW-l(^)) ^(^r"W(x),r"(^)-nx),...,x). ^ r(xo)). 

Proof. By Theorem 14.21 there exists a measure P on 9T(xo) such that 

ti(x) = / K{x,a) dO^a), (x G T(xo)). 
J9r(xo) 

Define the measure on Oaj^ by = P o <I). Then 



u{x) = / K(x,^>(xo,xi,...))(iz^(xo,xi,...). 

But i^(x,<I>(-)) is a multiple Cx of the characteristic function of Vj,„(x)(^x) r"('^)-^{x) x- From this, 
()4.13p follows immediately. □ 

Definition 4.6. A non-negative function on T(xo) is called additive if 

(4.14) u{x)= u{y), (xGr(xo)). 

r(y)=x 

Denote by 

Ty(")(x) = W{x)...W{r'''^{x)). 

Corollary 4.7. For any p-harmonic function u > there exists a unique additive function v such 
that 

(4.15) u(x) = —^-L_u{x), (x G r(xo)). 

Conversely, if u is given by ()4.15p /or some additive function v, then u is p-harmonic. 
Proof. The existence can be obtained from Theorem 14.51 by defining the additive function 

v{x) := iy{V^n(.)(^x),r"(-)-Hx),...,x)^ e T(xo)). 

Uniqueness is clear since W ^ on T{xq). 
For the converse, we compute 

V W{y)u{y) = V W{y) ^ nM-u 

,^x r^=x W{y)W{r{xi))...W{r^^y) i(y)) 
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" W{x)...W{r^i-)~\x)) "^^^ ^ 

r{y)=x 

□ 

Remark 4.8. Note that the function 2^0(2^) = W^"'^^^\x), x G T{xq) is an additive function. Indeed 

r(y)=x r{y)=x r{y)=x 

= Vr(x)...VF(r"(^)-i(x)) W{y) = iyo{x)-l. 

r(y)=x 

Therefore we have the fohowing corohary. 

Corollary 4.9. Every non-negative harmonic function is the quotient of two additive functions. 
Definition 4.10. A function U on T(xo) is cahed a QMF-weight if U > and 

Y U{y) = l, (xGr(xo)). 

r{y)=x 

Proposition 4.11. There exists a one-to-one correspondence between positive additive functions 
and positive QMF-weights on T(xo). 

For every positive additive function v on T{xq) the function 

U.{x) = ^^ ixeT{xo)\{xo}), 
v{r{x)) 

is a QMF-weight. 

Conversely, for every positive QMF-weight U , the function 

v{x) = u{xo)U{x)U{r{x)) . . . U{r<-^-\x)), {x G r(xo) \ {xq}), 

is additive, where i^(xo) is some fixed non-negative constant. 

Proof. We have 

r{y)=x r{y)=x r{y)=x 



The converse foUows as in Remark 14.81 □ 
Definition 4.12. Let be the set of points xq ^ X such that r"(2;o) is regular for all n > 0. 
Remark 4.13. Note that Xj. is invariant for both r and r~^. 



Proposition 4.14. Let h he an Rw-harmonic function on X^, i.e., equation ()2.6p is satisfied for 
X G Xr. Then for each xq G X^ there exists a unique additive function v^q on T{xq) such that 

(4.16) h{x) = "^^^ , {x G r{xo)). 

Moreover, the functions v^q are related by 

(4.17) Ur(^^^^){x) = W{xo)iy,^^,{x), (x G T(xo)). 
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Conversely, if Vxo o,n additive function on T{xq) for all xq € X^, and the functions satisfy the 
relation (j4.17p then the function h on defined by 

(4.18) h{x) = '^^^ = Mx), [x^Xr), 

is Rw -harmonic on X^. 

Proof. Since the restriction of h to T{xq) is p-harmonic, the existence and uniqueness of v^q such 
that ()4.16p is satisfied follows from Corollary 14.71 
We have for x £ T(xo), x is also in T(r(xo)) so 



Since nr(xp)(x) = n^Qix) + 1, and r""=o(^)(x) = xq, we have 

W{xo)iyxoix) = Mxo)ix)- 

For the converse, we compute 

r{y)=x r(y)=x r{y)=x 



□ 
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